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ABSTRACT 

We have conducted a deep (15 < r < 23), 20 night survey for transiting 
planets in the intermediate age (~ 550 Myr) open cluster M37 (NGC 2099) 
using the Megacam wide-field mosaic CCD camera on the 6.5m Multiple Mirror 
Telescope (MMT). In this paper we present a catalog and light curves for 1445 
variable stars; 1430 (99%) of these are new discoveries. We have discovered 20 
new eclipsing binaries and 31 new short period (P < 1 day) pulsating stars. 
The bulk of the variables are most likely rapidly rotating young low-mass stars, 
including a substantial number (> 500) that are members of the cluster. We 
identify and analyze five particularly interesting individual variables including a 
previously identified variable which we suggest is probably a hybrid 7 Doradus/5 
Scuti pulsator, two possible quiescent cataclysmic variables, a detached eclipsing 
binary (DEB) with at least one 7 Doradus pulsating component (only the second 
such variable found in an eclipsing binary), and a low mass (Mp ~ Ms ~ O.6M0) 
DEB that is a possible cluster member. A preliminary determination of the 
physical parameters for the DEB +7 Doradus system yields Mp = 1.58±O.O4M0, 
Ms = 1.58 ± O.O4M0, Rp = 1.39 ± O.OTi?© and Rs = 1.38 ± 0.07Rq. 
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Introduction 



This paper is the second in a series of papers on a deep survey for transiti ng planets in the 



open cluster M37 (NGC 2099) using the MMT telescope. In the first paper (IHartman et al. 



20071 . Paper I) we introduced the survey, described the spectroscopic and photometric obser- 
vations, determined the fundamental parameters (age, metallicity, distance and reddening) 
of the cluster, and obtained its mass and luminosity functions and radial density profile. As 
a result, we have a good understanding of the population of stars that we have observed. 
This paper focuses on the variable stars in our survey. We will discuss stellar rotation and 
transiting planets in future papers. 

As discussed in Paper I, there have been a number of surveys for transiting planets 
in galactic open clusters. The rationale behind conducting such a survey is that since the 
properties of stars in clusters are easier to determine en masse than they are for field stars, 
and since cluster stars have more uniform properties than field stars, one can determine in 
a relatively straightforward fashion from the survey the frequency of p lanets in the cluster, 
or at least place a meaningful upper li mit on it fe.g. iBurke et al. 1120061) . We were motivated 
in particular to conduct this survey by iPepper and Gaudil (l2005l . l2006l ) who suggested that 
it may be possible to find Neptune sized planets transiting solar-like stars by surveying an 
open cluster with a large telescope. These authors developed a formalism for estimating the 
expected yield from an open cluster transit survey as a function of the cluster's parameters 
and the survey strategy. Using this formalism we found that based on its location on the 
sky, distance, reddening, and richness, M37 is one of the clusters with the highest expected 
yield of Neptune-sized planets for a transit survey conducted with the MMT. 

While to date no confirmed transiting planet has been found in an open cluster, these 
surveys have produced a number of by-products unrelated to the study of planets. One 
of the primary by-products is the study of variable stars in the surveyed field. Studies 
of this sort have been published for the open clusters NGC 6791, NGC 2158, NGC 1245 



M34, NGC 2516, and NGC 2547 as a by-product of transit survevs flMocheiska et al.ll2002 



2004J : IPepper and Burke ll2006l : lDe Marchi et al.ll2007l : llrwin et al. Il2006l . l2007al ]bf). Similarly, 
variable stars have also been st udied in the globular clusters 47 Tucanae and u) Centauri a s 
a by-product of transit surveys (IWeldrake et al.ll2004l : IWeldrake. Sackett and Bridged 120071 ). 



The study of variable stars in stellar clusters is particularly interesting since one can 
put constraints on the age, metallicity, distance and reddening of stars that are cluster 



^Observations reported here were obtained at the MMT Observatory, a joint facihty of the Smithsonian 
Institution and the University of Arizona. 



- 3 - 



me mbers. This additional inform ation aUows for powerful tests of stellar pulsation models 
(e.g. 



[e.g. 



Frandsen and Arentoft 1119981) . or stellar evolution theory in the case of eclipsing binaries 



Southworth. Maxted. fc Smalley II2004I ). By identifying variables in differe nt clusters it 



Irwin et al- 



ls po s sible t o study the evolution of stellar propert ies including rotation (e.g 
20061 . l2007al jbl) and pulsation (e.g. iKrisciunas 1119981 ). Even if most stars in a given survey 
are not members of a cluster, identifying the variables is still interesting as it is a good tool 
for finding rare or previously unknown phenomena that deserve detailed study. 

While we were primarily concerned with maximizing the expected planet yield when 
choosing to survey M37, this cluster is also an interesting object to survey for variable stars 
both because it is rich for an open cluster (~ 4000 stars larger than 0.3Mq, Paper I) and 
because at an age of 550 Myr (Paper I) many of the stars are active enough to be detectable 
as variables with amplitudes greater than a few millimagnitudes. M37 has bee n the target 
of two previous surveys for variable sta r s (|Kiss et al. I l200ll : iKang et al. 1 120071 ) which have 
discovered 24 variables. The iKiss et al. I (120011 ) survey consisted of 658 i?-band observations 
taken over the course of ~ 1 month, and spanned a magnitude range of 12 < < 18, with 
a point-to-p o int pr ecision of 0.01 mag at the bright end and 0.03 mag at the faint end. The 
Kang et al. I (120071 ) survey consisted of 1220 V^-band observations taken over 12 nights in 
2004 and 3 nights in 2006, and was saturated for stars brighter than V = 13. Because our 
survey, consisting of ~ 4000 observations taken over ~ 1 month, and spanning a magnitude 
range of 15 < V < 23, with a point-to-point precision of ~ 0.002 mag at the bright end and 
0.05 mag at the faint end, has more observations with higher precision for fainter targets than 
either of these surveys, we have identified 1445 variables, of which 1430 are new discoveries. 

In the following section we will describe the observations and the reduction of the time- 
series data. In §3 we describe the selection of variable stars and present the catalog. In 
§4 we match to the previously identified variables. In §5 we discuss the global properties 
of the variables noting that a substantial fraction are likely to be rapidly rotating spotted 
stars. In §6 we analyze five particularly interesting variables including a probable hybrid 7 
Doradus/5 Scuti pulsator, two possible quiescent cataclysmic variables, a detached eclipsing 
binary with a 7 Doradus component and a low mass detached eclipsing binary with ~ O.6M0 
components that is a possible cluster member. Finally we summarize our results in §7. 
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Observations and Data Reduction 



2.1. Observations 



A detailed discussion of the observations was presented in Paper I, we provide only a 
brief overview here. The observations consist of both gri photometry for ~ 16000 stars, 
and r time-series pho tometry for ~ 23000 stars obtained with the Megacam mosaic imager 
jMcLeod et ^ boool ) on the 6.5 m MMT and high-resolution spectr oscopy of 127 stars 



obtai ned with the Hectochelle multi-fiber, high- dispersion spectrograph (jSzentgyorgyi et al. 
1998h on the MMT. 



The primary time-series photometric observations were done using the r filter and consist 
of ~ 4000 high quality images obtained over twenty four nights (including eight half nights) 
between December 21 ^ 2005 and Jan u ary 2 1, 2006. Where available we took BV photometry 



for the sources from iKalirai et al. I (120011). We also t ransfo rmed our ri photometry into 



Ic using the Ic photometry from iNilakshi and Sagar I (I20021) to defi ne the transformation. 
Finally we took Ks photometry from 2MASS (ISkrutskie et al. II2006I ) where available. 



The spectra were obtained on four separate nights between February 23, 2007 and 
March 12, 2007 and were used to measure Te//, [Fe/H], vsini and the radial velocity (RV) 
via cross-correlation. 



2.2. Image Subtraction Photometry 

The CCD calibration, point-spread-function fitting photometry and photometric cal- 
ibration routines were described in Paper I. To obtain light curves from the r' time se- 
ries observations of the pr imar y M37 field we used the image subtraction technique due to 
Alard and Lupton I (119981 ) and lAlard I (120001 ) . using a slightly modified version of Alard's 
Isis 2.1 package. Before performing photometry we construct a reference image for each 
chip from ~ 100 of the best seeing images. We choose one image as a template and register 
each image in the reference list to the template and then match the background and convolve 
it to the seeing and fiux scale of the template. The matched images are then mean-combined 
and processed through PSF fitting routines as described in Paper I. 



The reduction of each image in the time series involves the following steps: 



1. Remove cosmic rays using the extract program in Isis. 

2. Id entify stars in the image w ith the Source Extractor program (version 2.3.2) due 



to 



Bertin and Arnouts I (119961 ). The resulting star list is matched to the star list for 
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the reference image and a third order polynomial transformation from the reference 
coordinates to the image coordinates is determined. The reference image is then trans- 
formed to the image coordinate system using a cubic spline interpolation with the ISIS 
program interp. 

3. If the image has seeing better than 1" convolve the image with a Gaussian kernel to 
broaden the seeing up to 1". 

4. Adjust the FWHM values of the Gaussians used in the image subtraction kernel ac- 
cording to the value expected, assuming Gaussian PSFs, from the difference in seeing 
between the (broadened) image and the reference. We use the expected FWHM, half 
the expected FWHM and twice the expected FWHM. 

5. Use the mrj_phot program in IsiS to convolve the reference into the image, saving the 
residual image which is set to the flux scale of the reference. 

6. Transform the reference star list to the image coordinates and perform aperture pho- 
tometry on the residual image at the location of the transformed stars. To perform the 
aperture photometry we sum all pixels within a fixed aperture radius and divide by 
the sum of the reference PSF convolved with the kernel from mrj_phot over the same 
aperture radius. Here we use the PSF from the PSF fitting routines discussed in Paper 
I as the reference PSF to ensure proper flux scaling in converting the light curves from 
differential flux to magnitudes. 

The differential flux light curves obtained from the above pipeline are then converted to 
magnitudes using the reference PSF magnitude. We obtained differential light curves for a 
total of 23790 sources detected in the stacked r' images, of which 14750 arc detected in g, r 
and i, 7245 arc detected only in r and i, 141 are detected only in g and r and 1654 are only 
detected in r. 

The light curves are passed through two cleaning procedures to remove unreliable mea- 
surements. The first procedure is to remove 5cr outliers from the light curves, where a is 
the measured, rather than the expected, RMS of the light curve. This step is necessary to 
remove unidentified cosmic ray hits that fall within the aperture radius of a source, satellite 
trails or transient bad rows (which were seen to affect a number of images). All of these 
artifacts might affect only a few sources on a given image, but over the course of the run will 
affect almost every source at least a few times. Variable stars with sharp features, such as 
deep eclipses, may have reliable extremum measurements clipped by this procedure. There- 
fore, we investigate all identified variable star light curves on a case by case basis to recover 
reliable measurements that have been removed. 
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The second cleaning procedure is to remove points from images that are outhers in a 
large number of light curves. Letting x be the fraction of light curves for which a given 
image is a Sa outlier, we inspected the histogram of x and chose its cutoff value for each 
chip. The typical cutoff was between x = 5% — 10%. After both cleaning procedures, typical 
light curves contained ~ 4000 points. 

Figure [1] shows the Root Mean Square (RMS) of the light curves as a function of the r 
magnitude. The point-to-point precision of the light curves approaches 1 mmag at the bright 
end. This is approximately a factor of 5 poorer than one would expect based on the photon 
noise. Achieving better precision would require correcting for additional systematic varia- 
tions in the light curves. We will correct for these systematic variations before performing 
the search for transit candidates. 



3. Variable Star Selection 

We have used the following three different methods to identify periodic variable stars 
from our time series observations: 



1. A search for periodic, sinusoidal variations with th e Lomb-Scargle (L -S) algorithm 
Jlomb Ill976l : IScargle lEoS^ : IPress and Rvbicki lEoSol : IPress et al. Ill992h . 



2. A search for ge neral periodic variability usin g the analysis of var iance (AoV) peri- 
odogram due to ISchwarzenberg-Czerny I (Il989l ). as implemented by iDevor I (120051 ). 



3. A search for eclipses using the B ox-fitting Least Squares (BLS) algorithm due to 



Kovacs. Zucker and Mazeh I (120021 ). 



Note that the version of AoV that we implement here uses phase binning. ISchwarzenberg-Czerny 



( 119961 ) has also developed a version of AoV that uses triganometric orthogonal polynomials. 
This alternative version reduces to a form that is analogous to L-S when only one harmonic 
is used. Besides the difference in the variable used in the periodogram (AoV as opposed to 
power in L-S, note that there is a one-to-one correspondence between the two), there is a 
slight improvement in this version of AoV over L-S in that AoV also fits a constant term to 
the light curve rather than simply subtracting off the mean as is done in L-S. 

While there will be overlap between the results of the methods that we used, each 
algorithm tends to be sensitive to different types of variables. The more general tool, AoV 
with phase binning, should in principle identify variables found by the more restrictive L-S 
and BLS algorithms, but it may also be more sensitive to non-astrophysical variability due 
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to uncalibrated instrumental or atmospheric effects. As a result, AoV with phase binning 
cannot be used to search for sinusoidal or transit-like variability with amplitudes as low as 
that accessible to L-S and BLS. In the next three subsections we describe the results from 
using each selection tool, and in the last subsection we present the composite catalog of 
variables. 



3.1. L-S search 

We used the L-S algorithm to search for sinusoidal periodic variables with periods be- 
tween 0.01 and 20 days; the search was conducted with a sampling of 0.005 times the Nyquist 
sampling. Figure [2] shows the histogram of peak frequencies for all of our 23,790 light curves. 
Spikes in the histogram result from instrumental or atmospheric artifacts that systematically 
affect numerous light curves. Many of these occur at aliases of 1 sidereal day and result from 
artifacts such as the rotation of diffraction spikes from nearby stars or color dependent extinc- 
tion. Other artifacts, such as a series of observations where the image subtraction procedure 
did not work optimally (due, for example, to coma or other non-symmetric distortions of the 
PSF) may be phased at periods that are not aliases of a sidereal day. Since the vast majority 
of candidate variables identified at these frequencies are spurious detections, we reject any 
detection at one of these spikes. The rejection removes 8791 of the 23,790 light curves. The 
lower panel of figure [2] shows the histogram after removing these frequencies. The rise in the 
number of detections towards short frequencies results b oth from real variables with peri ods 



between 5 and 10 days and from ~ 1// red noise (e.g. iPont. Zucker and Queloz II2006I ) in 
the light curves. 

For each light curve we calculate the false alarm probability (FA P) of the per i od co rre- 



sponding to the peak in the periodogram using equation 13.8.7 from lPress et al. I (Il992l ): 

FAP(z) = 1 - (1 - e-^)*^, (1) 

where z is the peak value in the periodogram and M is the number of independent frequen- 
cies sampled, and which we approximate by the number of points times the ratio of the 
maximum frequency searched to the Nyquist frequency. For convenience we henceforth refer 
to log(-FAP) as FAP. The FAP provides a relative ordering of detection confidence among 
the light curves. How the formal FAP relates to the real FAP depends on the specific window 
function of the observations as well as on the noise distribution in the light curves. As the 
noise in our light curves is typically not frequency independent we do not expect the formal 
FAP to correspond to the real FAP, thus the selection threshhold in formal FAP is set by 
inspecting the light curves at a range of FAP values. 
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In figure [3] we plot the histogram of formal FAP values together with example light 
curves at a range of values. Only stars detected in all three filters (gri), with more than 
500 points in their light curves, and with r < 20 are included in the histogram. As seen in 
figure m we adopt a selection criterion of FAP < —150 for stars with P > 0.1 days and 
FAP < — 10 for stars with P < 0.1 days. We use a less severe selection threshold for short 
period stars since they are less affected by red noise. This selects 1533 candidate variables. 
After applying the selections described below in §3.5, 1369 of these candidates are included 
as variables in the final catalog. 

Note that while the formal false alarm probability for the long period detections is an 
astonishingly small 10~^^°, the actual false alarm probabilit y is much higher due to red noise . 



While metho ds such as the trend filtering algor ithm (TFA; iKovacs. Bakos and Noyes II2005I ) 



or SYSREM (ITamuz. Mazeh and Zucker II2005I ) may filter red noise from an ensemble of light 



curves, they are known to distort the signal of large amplitude variables and may reduce the 
sensitivity to these variables as well. Since we are primarily interested in large amplitude 
variables, we do not apply either of these filtering methods to our light curves. A filtering 
method will be applied to the light curves before searching for planetary transits and other 
low amplitude variables; the results from this will be presented in a separate paper. 



3.2. AoV search 

We used the AoV algorithm with phase binning to search for variables with periods 
between 0.1 and 20 days; a coarse scan was conducted at 0.01 times the Nyquist sampling 
with a fine tu ning search ar ound the highest peaks at 0.0005 times the Nyquist sampling. 



Following the iDevor I (120051 ) implementation, we normalize the AoV periodogram of each 
light curve to have zero mean and unit variance. As for L-S, we first reject spikes in the 
frequency histogram (removing 8968 light curves, figure [5]) and set the selection criterion by 
visual inspection of light curves. Figure [6] shows the histogram of normalized AoV values 
together with example light curves. We adopt a selection criterion of AoV > 3.75 which 
selects 982 candidate variables, 615 of these are included in the final catalog. 



3.3. BLS search 

We search for eclipsing variables using the BLS algorithm. The search is restricted to 
periods between 0.2 and 5.0 days as the detection efficiency drops for longer periods. As for 
L-S and AoV we first reject bad periods by searching for spikes in the histogram of detected 
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periods (5621 light curves are rejected; see figure [7]). 

Following iPont. Zucker and Queloz I J2006h we use as our selection statistic the signal 
to pink noise: 

= 2/ 1' 2/M (2) 

where 6 is the depth of the eclipse, nt is the number of points in the eclipse, Nt is the number 
of distinct eclipses sampled, o"^ is the white noise, and cr^ is the red noise at the time-scale 
of the eclipse. We subtract the best BLS model from each light curve and set equal to 
the standard deviation of the residual. We then bin the residual light curve in time with a 
binsize equal to the duration of the eclipse and set ar equal to its standard deviation. 

Figure [8] shows a histogram of the resulting SN values for all light curves, together with 
sample light curves at a range of SN values. We select 90 light curves with SN > 12.0 as 
candidate variables, 43 of these are included in the final catalog. 



3.4. Estimation of the Period Errors 

For the variables with periods detected by L-S or AoV, we measured the period un- 
certainty by bootstrap. For each bootstrap iteration we draw a sample of points from the 
light curve with replacement and then search for the period within a frequency of ±1/T 
from the period found for the original light curve, where T is the time-baseline of the light 
curve. The period uncertainty for a light curve is given by the standard deviation of 100 
bootstrap iterations. Figure [9] shows the period uncertainty of the variables as a function of 
their period. 



We can estimate the expected uncertainty on the period using equation 25 of lSchwarzenberg-Czerny 



tell): 



where P is the period, a is the noise in units of the signal amplitude (i.e., the inverse of 
the signal to noise), D is the average correlation time-scale for the residuals from the model 
periodic signal (equal to the sampling time for the case of pure gaussian white noise with a 
perfect sinusoidal signal), and T is the time base-line. This relation holds for T ^ P ^ D . 
Note that the period error does not include uncertainties due to aliasing - it describes how 
well you can centroid the peak in a periodogram after choosing amongst aliases. We choose 
by eye the best period for each light curve from among the possible aliases. 



We plot equation [3] on figure [9] for T = 31.1 days, assuming a typical signal to noise of 
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10 and a correlation time-scale of half an hour. Note that D is effectively a free parameter, 
since it is difficult to know a priori how well the AoV "model" will fit the actual light curves, 
and it is adjusted to provide a good fit to the short period errors. The errors from the 
bootstrap appear to be in reasonable agreement with theoretical expections, though only if 
we allow for time-correlated noise. The disagreement for the longer period detections is due 
to the violation of the assumption T ^ P. 



3.5. Combined Catalog of Variables 

The combined list of candidate variables contains 1963 sources. Of these, 38 were 
detected by all three methods, 543 were detected by AoV and L-S but not by BLS, 6 were 
detected by L-S and BLS but not by AoV, 16 were detected by AoV and BLS but not by 
L-S, 384 were detected only by AoV, 946 were detected only by L-S, and 30 were detected 
only by BLS. 

Before proceeding we must first correct for the blending of variability bet ween sources 



that results from performing aperture photometry on the subtracted images (e.g. lHartman et al 



20041 ) . To do this we identify groups of candidate variables that lie within 4" of each other (so 
that their PSFs may overlap significantly). We then choose from each group the source with 
the smallest average deviation of the centroid position measured on the subtracted images 
from the star's position measured on the reference image. After this procedure we are left 
with 1825 candidate variables. Requiring that the variables be detected in the g, r and i 
filters reduces this number to 1732. 

We then examine all candidate light curves by eye, choosing among the possible periods 
the one that appears to phase the light curve most cleanly, and rejecting false positive 
detections. In addition to inspecting the hght curves we also inspected the position of each 
star on the reference mosaic image rejecting stars on bad columns or in the wings of very 
bright stars. After this procedure we are left with 1445 variables, of which 1380 have periods 
(note that variables without periods correspond to light curves which showed believable 
variability but for which the period returned by the selection tools did not appear to be 
correct). In figure [10] we show examples of some of the 287 light curves which were rejected 
by eye. Note that some of these may well be real variables, and unfortunately the selection 
process is fairly subjective, however it is very difficult to develop an automatic selection 
routine that is completely robust against false positives, but does not reject a significant 
number of true variables, when dealing with real data. 



Table [T] gives the first few entries in our catalog of variable stars, the full catalog is 
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available in the e lectronic edition of the journal. We c ontinue the number ing system from 
Kiss et al. I koQlh . relabeling KVl through KV17 from iKang et al. I J2007h as V8-V24 and 
labeling new variables begining with V25. The new variables are sorted by RA. For com- 
pleteness we include unmatched known variables, though we caution that the coordinates for 
V8-10, V15-16, V18, and V23 are uncertain since we did not observe these stars and there 



appears to be an error in the table listing the variable star coordinates in lKang et al. I (120071 ) 
(see next subsection). Light curves for all variable stars are available on request. Table [2] 
provides a summary of the variable star catalog. 

We also inspected the light curves of the 94 variable candidates that were not detected 
in all three filters, we selected 29 of these as real variables based on their light curves. Table [3] 
lists the positions, magnitudes and periods for these suspected variables. 



Match to Known Variables 



There are a total of twenty four varia bles that have bee n previously identified in the 
field of M37. Of these, seve n were found by iKiss et al. I (1200 ll ) and an additional seventeen 
were found by iKang et al. I (120071 ). Thirteen of these variables are eclipsing binaries, ten 
show 5-Scuti type variations and one is an RRc variable. 

We matched the previous variable star catalogs to the list of 23,790 sources that we 
obtained light curves for using an astrometric matching radius of 0'.'5. Fifteen of the known 
variables matched to sources in our input photometry list. Six of the unmatched variables 
(Vl-2, V8-10 and V18) are saturated in our observations, the other three unmatched sources 
(VI 5- 16, and V23) fell on chip gaps and were not observed. We independently select all 
fifteen matched sources as variables. In matching to the known variables we found that 
there appears to be a discrepan c y betw een the coordinates for several of the KV variables 
listed in table 1 of iKang et al. I (120071 ) and the light curves presented in that paper. We 
came to this conclusion by comparing the light curves and periods of the matched sources 
with the light curves in the paper. We revise the coordinates in our own catalog following 
the naming convention used for the light curves in iKang et al. I ( 120071 ). 



Figure [TT] shows our light curves for the fifte en variables that we recover together with 
the published photom etry from iKiss et al. I (120011 ). For the variables V3-V7 we combine our 
photometry with the iKiss et al. I (120011 ) photometry which was taken between Dec. 1999, 
and Feb. 2000 to improve the precision of the periods. We will now discuss a few of these 
variables in detail, though we leave the very interesting star V17 for §6. 
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4.1. V3 and V4 



The contact e clipsing binaries V S and V4 both exhibit a pronounced O'Connell effect 
(unequal maxima; lO'Connell I Il95ll ). The origin of this effect is still uncertain, though 
a number of theo ries have been put forward to explain it including stellar spot activity 
( Binnendiik I Il960l ) . the formation of a hot spot on one of the stars due to mass exchange 
( Shaw 1994[). or an unequal in teraction between the two components and any circumstellar 
material (ILiu and Yang II2003I ). In any of these theories, one would expect the effect to be 
variable. V4 does show some evidence for variation of ~ 0.01 mag between the 2000 and 
2006 epochs as well as within the 2006 epoch as seen in figure [TTl V3 on the other hand, 
does not show evidence of variation. Note that four of the newly discovered contact binaries 
also exhibit the O'Connell effect (V37, V706, V936, and V1194, see figure US]). 



4.2. V19 



Kang et al. I (120071 ) classify this variable as a possible W UMa type eclipsing binary. 



though the classification is unsecure due to a lack of phase coverage. We find that the light 
curve shows substantial amplitude variations which make it most likely a fundamental mode 
RR Lyrae pulsator with the Blazhko effect. 



5. Properties of the Variables 

In figure [T2] we show the location of all the variables identified by this survey on gr and gi 
CMDs. It is immediately apparent that the variables lie preferentially along the cluster main 
sequence. Most of these are rapidly rotating stars (relative to the Sun) with surface brightness 
inhomogeneities that are detectable at our level of photometric precision. We will discuss 
these variables in detail in a separate paper (Hartman et al. 2007, in preparation). Given 
the relatively young age of the cluster (550 Myr), the fact that we detect these variations 
for so many cluster stars is not surprising. 

From inspecting the light curves of all the variables identified by our survey, we classify 
27 of the variables as eclipsing binaries (or ellipsoidal variables) and 39 of the variables 
as pulsating variables with periods less than 1 day. Including the 4 previously identified 
eclipsing binaries and 5 previously identified pulsating stars that we have not observed, 
there are a total of 31 eclipsing binaries and 44 pulsating stars known in the field of this 
cluster. We show the locations of these variables on the gr and gi CMDs in figure [13] and 
show the phased light curves for the newly discovered variables in figures [T3] and [TS] As the 
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red-edge of the (5-Scuti insta bility strip intersects the M37 main sequence at \^ ~ 14.4 or 



r ~ 14.2 (iKang et ah 1120071 ). we expect that any 5-Scuti variables within this cluster will 
be saturated in our time series observations. The one or two pulsating variables in figure [12] 
that lie near the cluster main sequence are thus likely to be background disk stars. There 
are six eclipsing binaries that lie near the cluster main sequence, these include V3, V4, V692, 
V855, V1028 and V1380. The binaries V3, V4 and V855 are contact systems, while V692, 
V1028 and V1380 are detached systems. We will discuss V1028 and V1380 in greater detail 
in §6.4 and §6.5. 

We have obtained spectra for the system V692 and detected it as a single lined spectro- 
scopic binary. The cross-correlation routine gives a temperature for this system of ~ 7000 
K which is too hot for it to be a member of the cluster given its position on the CMD, but 
consistent with it being a member of the background field. Based on the radial velocity 
curve, we estimate that the systemic RV is ~ —2 km/s compared with 9.4 ± 1.4 km/s for 
the cluster (Paper I). We conclude that this system is not a member of the cluster, and will 
not investigate it any further. 

Figure [11] shows the period-r relation for the variables. In two of the panels we split the 
variables into two groups, those which match to photometrically selected cluster members 
(see Paper I), and those which lie well away from the cluster main sequence. There is a strong 
correlation between period and magnitude for the cluster members with a break into two 
groups at r ~ 18. We are most likely seeing the relation between rotation period and stellar 
mass for the y oung cluster stars. T he relation is quite similar to that seen, for example, in 



the Hyades by iRadick et al. I (119871 ). and will be discussed in greater detail elsewhere. The 
non-cluster members appear to show no clear correlation between period and magnitude. 
This is not unexpected, though, since the field stars lie over a large range of distances from 
the Earth. 

We will now discuss a few of the populations of field variables. 



5.1. Short Period Pulsators 



The majority of the short-period pulsating variables appear to lie in the range r < 19, 
0.2 < g — r < 1. These stars have pulsation periods on the order of 0.1 days and have 
0.45 < B — V < 1.18, they are thus likely to be background 5-Scuti or SX Phe (Population 
II 5-Scuti) stars largely seen through E{B — V) > 0.4. 

Fundamental mode (5-Scuti variables exhibit a well -known period- lura i nosity -color (PLC) 
relation that we can use to determine their distances. iMcNamara et al. I (120001 ) give the fol- 
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lowing period-color and period-luminosity relations: 

{V-Ic)o = 0.195 logP + 0.514, 
My = -3.725 log P- 1.933. 



(4) 
(5) 



The first of these equations can be used to determine the E{V — Ic) reddening, then assuming 
E{y — Ic) = 1.5QE{B — V) (Paper I) and Ay = 3.1E{B — V) we can determine the distance 
using equation O 

We identify four stars, based on their periods, amplitudes and light curve shapes, that 
are quite likely fundamental mode pulsators. We give the identificatio ns, colors, periods and 
inferred distances in table HI t he BV photometry f or V6 comes from iKang et al. I (120071 ) as 
this star was not observed b v iKalirai et al. I (1200 ll ). We include V5 in this list, though we 
note that iKang et al. I (120071 ) have argued that this star may be an RRc variable rather than 
a fundamental mode 5-Scuti. Figure [T7| shows the distance and extinction for these four 
stars as well as for a number of eclipsing binaries discussed below. One star lies well beyond 
the disk of the galaxy and is likely to be a halo SX Phe star. The total extinction in this 

y larger than the 
) extinction maps 



direction appears to be Ay ~ 2 ma g, or E(B — V) = 0.65 which is slight 
value of E{B — V) = 0.56 from the ISchlegel. Finkbeiner and Davis I (|l998 



(note that the galactic latitude of the field is 
(119981 ) map is unreliable here). 



3°, so the 



Schlegel. Finkbeiner and Davis 



5.2. Eclipsing Binaries 

As seen in figure [T5|, many of the eclipsing binaries are contact, W UMa type systems. 
Similar to (5-Scuti variables, these binaries exhibit a PLC relation that we can use to de- 



termin e their distances. We adopt the following PLC relation from iRucinski and Duerbeck 



(119971 ): 

My = -2.38 log P + 4.26(5 -1/)o + 0.28 (6) 

where we assume solar metallicity, and note that the resulting uncertainty on the distance 
modulus for each source is ~ 0.5 mag. 

The reddening of each system can be estimated from the B — V-V — Ic color-color 
plot assuming E{V — Ic) = l.56E{B — V) and using the fiducial relation for the cluster. 
Assuming Ay = 3.1E{B — V), we can then find the distance to the system from equation El 
Table [S] gives the identifications, colors, periods and inferred distances and extinctions for 
10 candidate contact binary systems. We do not include two systems which appear to have 
unreliable photometry (one has a negative reddening, the other yields Ay > 4.0 mag), and 
one star which has dereddened colors that are blueward of the cluster turnoff. 
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Figure [T7I shows Ay as a function of distance for the 10 systems hsted in table [51 The W 
UMa systems appear to show more scatter than the 5-Scuti stars which may be attributed 
in part to the 0.5 mag uncertainty in the P — L relation for these stars and in part to the 
greater uncertainty in their colors due to their larger amplitude variability. There is also 
a hint that the W UMa systems have systematically greater extinction at a given distance 
than what we infer from the (5-Scuti stars. It is unclear what the source of this systematic 
discrepancy might be. 



We also show the extinction vs. distance from the Besangon model ( iRobin et al. II2003I ) 
for the galactic latitude/longitude of the field {I = 177.6°, b = 3.1°), assuming the default 
interstellar extinction of 0.7 mag/kpc. While the value of 0.7 mag/kpc appears to reproduce 
the extinction of the cluster, the model extinction/distance relation cuts off sharply at the 
scale-height of the galactic disc (~ 140 pc) whereas the observed relation cuts off at about 
twice that distance. The model assumes that the dust follows a simple exponential density 
profile in both the radial and vertical directions; discrete clouds along the line of sight with 
higher than average density may cause the total extinction in this field to be greater than 
expected. 



5.3. Rotating Field Stars 

There are ~ 500 variable stars that do not appear to be members of the cluster, are 
not eclipsing binaries, and are not 5-Scuti or RR Lyrae type pulsators. We argue that most 
of these stars are likely to be rotating variable main sequence stars. Figure [18] shows light 
curves for 16 random variables of this type. 

In figure [T9] we plot period against V — Ic for the 481 of these variables that have period 
determinations and V photometry. For stars with 1 < V — Ic < 2.3 the periods appear 
to scatter between 0.1 and 20 days, with the highest density between 3 and 12 days, while 
stars with V — Ic < 1 or V — Ic > 2.3 tend to have shorter periods (< 2 days for stars with 
V — Ic > 2.3). The trend toward shorter periods among the bluer stars may either indicate 
a population of pulsators (e.g. 7 Doradus stars) or it may indicate shorter rotation periods 
for the hotter stars. From a visual inspection of the light curves alone it is impossible to 
distinguish between these two cases. For the redder stars, however, the trend toward shorter 
periods is probably due to the onset of short period rotation for cool stars. Note that for 
the cluster, V — Ic = 2.3 corresponds roughly to r ~ 20, which, as seen in the center 
panel of figure [16], is where the stars appear to transition to a population of predominately 
short period rotators. The scatter between 3 and 12 days is also consistent with the range of 
rotation periods seen for brighter stars, keeping in mind that as the rotation periods increase 
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the amplitudes decrease (thus making longer period rotators harder to detect) and that the 
rotation periods generally evolve over time from rapid rotati on in the pre-main sequence 



stage to month long rotations for stars the age of the sun (e.g. Ilrwin et al. Il2007al ). 



One way to test the consistency of our hypothesis that these variables are rapidly rotat- 
ing, heavily spotted stars is to compare the amplitudes of their light curves to the Rossby 
number [Ro] the ratio of the rotation period to the convective turnover time-scale). Empir- 
ically there is a st rong anti-correlati on between the Rossby number and the x-ray activity 
of low-mass stars (INoyes et al.lll984l ). and a similar anti-corr elation appears to hold for the 



light curve amphtudes of spotted rotating stars as well (e.g. 
200lh . 



Messina. Rodono and Guinan 



The turnover time-scale can be estimated from the empirical relation given by lNoyes et al. 



(119841 ). To apply this relation we must first de-redden our stars. As for the W UMa stars, 
this can be done using the B — V vs. V — Ic color-color plot together with the fiducial 
de-reddened relation for the cluster (figure [20]) . Because the reddening vector for stars with 
B — V > 1.2 01 V — Ic > 1.8 may intersect the fiducial curve at more than one location, 
we select stars bluer than these limits. We also adopt a cut on the inferred E{B — V) of 
< E{B — V)< 0.7 and a magnitude cut of 15 < V < 22. To avoid including stars without 
surface convection zones we also only consider stars with {B — V)q > 0.5. This last cut 
rejects 27 stars, many of which may be 7 Doradus variables. 



In figure ED we show the relation between Rg and the peak-to-peak r-band amplitude 
(Ar) for the 232 stars which pass the above cuts (we also reject one star with ~ 2 as it 
is an outlier from the rest of the observed stars). By eye there is a clear anti-correlation. 
To show that this is not due to a selection effect, we also plot for each star the estimated 
minimum amplitude that the light curve could have had and still have passed the variable 
star selection criteria. The methods for calculating the amplitude of the light curves and 
estimating the minimum detectable amplitudes will be described elsewhere (Paper 111). To 
evaluate the signficance of the apparent correlation we calculate Kendell's normalized rank- 
based r_corre]ation_co the case of data suffering a one-sided truncation 
( lTsailll990l : lEfron and Petrosianlll992l . ll999l ). We find r = —2316 meaning the null hypothesis 
of no correlation is formally rejected at the 2316cr level. This correlation is stronger than what 
we find for the period- (r = —1728) or for the {B — V)Q-Ar (r = 565) correlations. The 
anti-correlation between Ro and Ar is consistent with the hypothesis that rotation coupled 
with surface brightness inhomogeneities is the cause of the observed brightness variations for 
these stars. 
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6. Individual Variable Stars 

So far we have focused primarily on groups of variables. In this subsection we will 
investigate a few particularly interesting individual variables, including a hybrid 7 Doradus/(5 
Scuti pulsator (V17), a blue periodic variable star with g — r < 0. (V1148), a blue flaring 
variable (V898), and two interesting detached eclipsing binaries (V1380 and V1028). 



6.1. V17 



Kang et al. I (120071 ) emphasize the curious nature of the variable VI 7. They find two 
independent frequencies in its light curve: one at /i = 22.796 day~^ with a semi-amplitude 
of 32.9 mmag in V, and a sec ond much lower fre quency of /2 = 2.396 day~^ with a semi- 
amplitude of 24.2 mmag in V. iKang et al. I (120071 ) give two possible interpretations for this 
variable, that it is a hybrid 7 Doradus/5 Scuti type pulsator, or that it is a close semi- 
detached binary system with a 6 Scuti primary and a late-type giant secondary showing 
ellipsoidal variations. They claim that the 2MASS color for the object suggests a late G or 
early K spectral type leading them to conclude that the binary hypothesis is more likely. 

Note that 7 Doradus pulsators are early F-type main sequence stars that have periods 
in the range of 0.4-3 days, anaplitudes of a few percent and typically oscillate in multiple 
frequencies (IKrisciunas Ill998l : iKaye et al. Ill999l ). The variations in these stars are thought 
to be due to high-order g-mode pulsations. The instability strip for this type of variable 
extends from 7200 — 7550 K on the zero -age main sequence (ZAMS) and from 6900 — 7400 
near the end of the main sequence phase (IDupret et al. 1120041 ). S Scuti pulsators, on the other 
hand, are late A or early F-type stars that lie in the intersection of the classical instability 
strip and the main sequence. These stars have periods less than 0.3 days and amplitudes 
ranging from a few mmag to almost a full magnitude. The variations in this case are due 
to p-mode pulsations. The 6 Scuti and 7 Doradus instability strips overlap so that hybrid 
stars, exhibiting both modes of pulsation, may exist. 

To aid in understanding this interesting object we have conducted a multi-frequency 
analysis of our r light c urve using the Discrete Fourier Tran sform (DFT) with the CLEAN de- 
convolution algorithm (IRoberts. Lehar and Dreher 1119871 ). Figure [22] shows the DFT power 
spectrum before and after implementing CLEAN; we also show the power spectrum of the 
window function for reference. We confirm the high-frequency signal at /i = 22.796 day~^ 
and we find four additional low frequencies of /2 = 1.418 day~^, /s = 2.263 day~^, = 2.054 
day~^, and /s = 0.823 day~^. None of these frequencies are harmonics of each other. The 
error on all of these measurements is ~ 0.003 day~^. 
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After identifying the frequencies using DFT/CLEAN we fit a series of cosinusoids to the 
hght curve of the form: 



+ ^ tti cos(27r/i(t - to) + "^i 



(7) 



to get the semi-amphtude and phase of each mode. For consistency with iKang et al. I (120071 ) 
we adopt to = 2453000.0. We find semi-amphtudes of ai = 26.4 mmag, 02 = 13.7 mmag, 
03 = 12.8 mmag, 04 = 12.9 mmag, and 05 = 11.2 mmag, and phases $1 = 2.854, $2 = 1.132, 
$3 = 4.866, $4 = 4.757 , $5 = 5.477. Alternatively, if we assume only the two frequencies 
from lKang et al. I (120071 ) of fi = 22.796 day"^ and /2 = 2.396 day~^ we find ai = 27.2 mmag, 
a2 = 8.0 mmag, $1 = 2.144, and $2 = 1.996. 

The presence of multiple frequencies in the range 0.8 < / < 2.3 indicates that the 
low frequency variations are not due to ellipsoidal variability, and is instead suggestive of 7 
Doradus pulsations. This is not in confl i ct wit h the color of the star since the late-G/early-K 
spect ral type estirnate by iKang et al. I (120071 ) assumes no reddening. Taking B — V = 0.6 
from iKalirai et al. I (120011 ) and noting that the reddening to the source must be at least the 
cluster value of E{B — V) > 0.23 given that the star lies below the cluster main sequence 
on the CMD, we flnd {B — V)d < 0.37 which, for a dwarf star, corresponds to a spectral 



type earlier than F5 (e.g. ISchmidt-Kaler Ill982l ). We conclude therefore that the star may 



very well be a hybrid 7 Doradus/(5 Scuti pulsator. These stars are pa rticularly interesti ng 
for astroseismology and to date only a handful of them are known (see iKing et al. 1 12006| ) . 



6.2. The Hot Variable V1148 

As seen from its position on the CMD, star V1148 is an outlier with respect to the general 
population of variables found in our survey (flgure [T3|) . The star has g — r = —0.396 ±0.013, 
r-i = -0.247 ± 0.015, r = 19.97 ± 0.012, B-V = -0.071 ± 0.004 and V = 19.770 ± 0.003. 
This places the star roughly 2 magnitudes above the cluster white dwarfs. As seen in flgure [H] 
the star shows saw-tooth like variability with a period of 0.157728 ± 0.000010 days and an 
amplitude of 0.07 mag in r. This variable is likely to be a quiescent cataclysmic variable, 
though we will consider several hypotheses below. 

While the position of the star on the CMD is consistent with it being a foreground 
white dwarf, the period is nearly an order of magnitude longer than the longest seen for 
white dwarf pulsations (the longest periods being 5-50 minutes for GW Vir stars whic h 



are still embedded in planetary nebulae, see ICorsico. Althaus and Miller Bertolami II2006I ). 



Such a long period could be possible for a binary system, with the variations resulting from 
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the reflection and/or ellipticity effects, or from an accretion hot-spot on the white dwarf. 
The sawtooth appearance of the light curve would be most consistent with the accretion 
hypothesis, though we note that the accretion rate would be quite stable over the course of 
a month if this hypothesis is correct. 



Since the star cannot have {B — V)o < —0.33 (jSchmidt-Kaler I Il982l ). the maximum 
reddening that it could be seen through is E{B — V) = 0.26 {Ay = 0.8). From figure \T7\ 
we can estimate that the star must be closer than ~ 5 kpc assuming that the star lies 
along an average line of sight. We can therefore put a limit on its absolute magnitude of 
My > 5.5. This rules out a main sequence star which must have My < 1.1 given the 



color (jSchmidt-Kaler 1119821 ). Even if there were no extinction, if the star were on the main 
sequence it would have to be at least ~ 50 kpc away, placing it well beyond the disk of the 
galaxy, or ~ 300 kpc assuming at least the reddening of the cluster. It is possible that the 
object could be a hypervelocity B star wh ich has been ejecte d from the galaxy, these have 
been found at distances of 50-100 kpc (e.g. iBrown et al.ll2007l ): though if this is the case the 
extinction to the star would have to be anomalously low. It is interesting to note, though, 
that the light curve shape, period, a mplitude and color are all consis tent with the star being 
a main sequence /3-Cephei variable (jSterken and Jerzykiewicz Ill993l ). 



It is possible that the source is a sub-dw arf B (sdB) star, though the typical a bsolute 
magnitude of an sdB is My ~ 4.6 or brighter (IMoehler. Heber and Rupprecht 1119971 ) which 
is slightly brighter than the limit we set based on the maximum extinction for the source. 
There is a class of puls ating sdB stars, ca lled PG 1716 stars, with periods comparable to 
what we see for VI 148 (iGreen et al. 1120031 ). These stars pulsate in high-order g-modes and 
typically have very low amplitudes (1 mmag or less), with the highest observed amplitudes 
of ~ 4 mmag (IKilkenny et al. 1120061 ). If V1148 is a member of this class, then its amplitude 
is more than an order of magnitude larger than the other members of its class. Such a large 
amplitude would be unusual for a high order g-mode pulsator since it would have to be 
pulsating in a very low degree / mode. Alte rnatively the star might be a short-period sdB 



in a binary system (e.g. lO 'Toole et al.ll2006l ). 



Alternatively, the star might be a blend between a white dwarf or an sdB star and a 
much fainter variable. If the variable is a background 5-Scuti or SX Phe star, the period 
and light curve shape suggest that it would be a fundamental mode pulsator. In that case 
we can estimate its intrinsic color and magnitude using equations H] and [51 This yields 
{V — Ic)o = 0.36, and My = 1.05, or B — V ^ 0.3. Assuming a minimum apparent B — V 
of —0.33 for the white dwarf/sdB star, and an extinction of Ay = 2.0 mag to the (5-Scuti, 
the (5-Scuti must have ^ > 21, and be d > 40 kpc away. As seen in figure [17] there does 
appear to be an unblended SX Phe star at d ~ 33 kpc, and it is possible that there could be 
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an even more distant one that happens to he along the hne of sight to a white dwarf. 

The chances for a random blend of this nature though are rather low. Out of the 16431 
point sources detected in our field with r and i measurements and with r < 24, there are 
128 pairs of objects that are separated by 0.1 — 1". We do not detect any evidence that the 
object is anything other than a point source, so if the objects are within two magnitudes 
of each other then we estimate that they must be separated by O'.'l at the most. We would 
expect ~ 1 such random pairing out of all point sources in our field. The chances that one of 
the paired objects is a white dwarf and the other is a distant halo SX Phe star (both fairly 
rare in our field) is thus very low. 

Finally, we note that the variable star may be a short period rotating K or M dwarf 
(§5.3) that is in a binary system with a white dwarf or sdB star. Such a pairing is more 
likely than the chance alignment of a white dwarf with a distant halo SX Phe star. A period 
of 0.157 days would be on the short side for such a star (see figure [19]), though double this 
period would also be permissible since the star may have multiple spots. The shape of the 
light curve is similar to that seen for other variables in the cluster. The spotted cluster stars 
have amplitudes Ar < 0.2 mag; taking this as the maximum amplitude for a spotted star 
in the r-band we put an upper limit on the difference in magnitude between the companion 
and the white dwarf of rr. nmp — Twn < 0-6 ma g. The minimum g — r value that a star can 



have is g — r = —0.6 (see iGirardi et al.ll2004l ). so to reproduce the color of the system the 



companion would have to have g — r < 0.09 (i.e. FO or earlier). Therefore, this hypothesis 
appears to be inconsistent with the observations. 

The most likely explanation for this variable is that it is foreground white dwarf with 
a close M dwarf companion, with the variations due to an accretion hot-spot on the surface 
of the white dwarf. Understanding the nature of this variable will require spectroscopy to 
determine log{g) and T^ff of the hot star and to check for evidence of accretion and/or a 
blend with a cooler star. 



6.3. V898 

The light curve of the variable V898 (figure [23]) is quite unique relative to the other 
variables in the survey. The light curve shows continuous non-periodic variability with a 
number of relatively small amplitude outbursts (0.1 — 0.2 mag). The typical duration of these 
outbursts is 0.5 days with a time-scale between outbursts of 1-2 days. The variation seen i n 



this source is similar to the flickering seen in many cataclysmic variables (see lDuerbeck Ill996l ) 



confirming this identification would require additional observations including spectroscopy. 
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The source matches to 2mass05523126+3236324 and has B - V = 0.604, V = 16.090, 
g-r = 0.474, r-i = 0.208, r = 15.876, J-Ks = 0.311, H - Ks = 0.110 and fsTs = 14.374. 
We estimate that the reddening to the object is E{B — ~ 0.29. 



6.4. Detached Eclipsing Binary V1380 

The echpsing binary V1380 is one of three detached echpsing binaries that he near 
the cluster main sequence. V1380 matches to 2mass05531108+3224434, it has an orbital 
period of 2.1916 ± 0.0016 days, and photometry B ~V = 0.594, V = 14.952, g-r = 0.42, 
r-i = 0.15, r = 14.79, J - Ks = 0.331, H - Ks = 0.084, and Ks = 13.331. The source is 
just at our saturation threshhold, and in many of the images it is saturated, so there are only 
1111 points in its light curve (a factor of ~ 4 fewer than most light curves). If the system is 
a member of the cluster, its primary star would have M ~ I.2M0 based on its r magnitude 
and assuming equal mass components. We note, however, that the star does appear to be 
slightly below the main sequence, so its photometric membership is questionable. 

To determine the period of the out of eclipse variations, we remove the in-eclipse points 
from the light curve and, using L-S, find a period of 0.939 ± 0.01 days with a peak-to-peak 
amplitude of 0.019 mag. We subtract a sinusoid fit to this signal and find in the residual 
light curve an additional period of 1.166 ± 0.01 days with an amplitude of 0.006 mag. To 
verify the detection of both periods, we have also used DFT/CLEAN to search for multiple 
periods. Figure [2l] shows the window power spectrum, the raw DFT power spectrum of the 
light curve, and the power spectrum after using CLEAN. Using DFT/CLEAN we find two 
strong period detections at 0.941 ± 0.011 days and 1.174 ± 0.013 days. In figure [231 we show 
the original light curve phased at the orbital period of the system, the out of eclipse points 
phased at 0.941 days, the residual light curve after subtracting the above period phased at 
1.174 days, and the full light curve after subtracting a double sinusoid fit to the out of eclipse 
points. The model is subtracted through the eclipses, though this is technically incorrect. 
We adopt the DFT/CLEAN periods for the rest of the discussion. 

The additional periods are most likely not the rotation periods of the two binary com- 
ponents. The time-scale for the synchronization of rotation periods wit h the orbita l period 



is much shorter than the time-scale for the circularization of the orbit (jZahn 1119771 ). Since 
the orbit appears to be circular (see below), we expect that it should also be synchronized. 
We will show that the additional periods are most likely due to one or both of the stars 
being a 7-Doradus type pulsator. 



We will now determine the physical masses and radii of the components of this binary 
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system. To do this we will combine radial velocity curves with our light curve. 

This binary was among the sources that we obtained spectra for with Hectochelle. The 
spectroscopic classification procedure discussed in Paper I yields a temperature of T^ff = 
6600±400 K, and v sini = 21.5±0.7 km/s, these values are an average of the two stars. The 
metallicity is unconstrained so we assume solar metallicity. Only the spectra for the first two 
nights are used in the classification since the two components were not well separated in the 
last two nights; classifying spectra from the last two nights yields an artificially high value 
for V sin i. For the cluster reddening of E{B — V)= 0.227 (Paper I), the expected color of the 
system is B — V = 0.65 ± 0.1 which is consistent with the observed value of B — V = 0.594. 

We combine the spectra from each of the four usable nights and calculate the cross- 
correlation as a function of radial velocity (RV) with xcsao using the best template from 
the classification procedure. The cross-correlation plots show two clear peaks, indicating 
that we do detect the system as a double-lined s pectroscopic binary. We th en use the TwO- 



Dimensional CORrelation algorithm (TODCOR; IZucker and Mazeh Ill994j ) to determine the 
velocities for each component. We use the same template for both stars as the nearly equal 
depth eclipses means that the stars have similar surface temperatures. Table lists the 
radial velocity and light ratio measurements from the spectra. The average light ratio in the 
RV31 filter is L2/L1 = 0.992 ± 0.025. 

To combine the radial velocity curves with the light curve we need to improve the period 
that is derived from the light curve. To do this we fit a zero-eccentricity orbital solution to 
the radial velocity curves, assuming the period determined from the light curve, and allowing 
for a shift in phase. From the phase shift we calculate a single "observed" primary eclipse 
time at the epoch of the RV observations. We then determine eclipse minimu m times frorn 



the light curve using the Detached Eclipsing Binary Light curve fitter (DEBiL: iDevor II2005I ) 
and fit for the ephemeris using the equation: 

T^in = HJDo + nP. (8) 

Here T^jn is the observed minimum time, HJDq is the reference primary eclipse epoch, P 
is the orbital period and n is the eclipse number. We find HJDq = 2453724.5306 ± 0.0004 
and P = 2.19258 ± 0.00004 days, where the errors are the standard errors from the the least 
squares fit. The observed eclipse timing measurements are listed in table [3 We include 
secondary eclipses in the fit and assume zero eccentricity, which we justify below. 

Using the above ephemeris, we fit a zero-eccentricity orbit to the RV curves finding Ki = 
118. 3± 1.3 km/s, K2 = 118. 3± 1.3 km/s and 7 = 5.3±0.7 km/s. The mass ratio is thus q = 
l.OOiO.Ol, so the two components have equal mass to within the precision of our observa tions. 



This is another example of a twin binary star system (jPinsonneault and Stanek II2006I ). The 
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systemic radial velocity 7 is more than 3cr below the mean cluster value of 9.4 ± 1.2 km/s 
(Paper I). Figure [26] shows the radial velocity observations with the best fit orbit. 



S ince the system is well detached, we model the light curve usin g the JKTEBOP pro- 



gram (ISouthworth. Maxted. fc Smallev 112004: ISouthworth et al. 



Eclip sing Binary Orbit Program (EBOP: lPopper and Etzel Ill981 



2004 ) which is based on the 



Etzel Ill981: Nelson and Davis 



19721 ). but includes more sophisticated minimization and error analysis routines. We assume 
a linear limb darkening law, adopting u = 0.564 for both the primary and secondary stars. 
This coefficie nt is appropr iate for the the r filter with Tg// = 6600 K, log{g) = 4.5 and 
[M/H] = 0.0 (IClaret II2004I ). We use a gravity darkening exponent of 1.0, which is the value 
expected for a radiative envelope star, though the effect is negligible in this system. The 
resulting parameters are given in table El where the l isted errors are the la errors from 
100 Monte Carlo simulations. ISouthworth et al. I (120051 ) have shown that the Monte Carlo 
routine in JKTEBOP provides robust error estimates, though we note that uncertainties in 
the limb darkening and reflection laws are not included in the Monte Carlo simulations. In 
conducting the simulations we set the error on each point in the light curve equal to 3.9 
mmag, which is the RMS of the residual from the best flt. The best flt model to the light 
curve is shown in flgure [271 

The ratio of the radii, k, is poorly constrained from the light curve. A better constraint 
could be obtained by including the spectroscopic constraint on the luminosity ratio, however 
since the spectra were obtained over a signiflcantly different wavelength range from the light 
curve, including this constraint is not straightforward and beyond the scope of this paper. 

To put a limit on the eccentricity of the system, we have also flt the light curve allowing 
ecosa; and esina; to vary. We flnd e = 0.003±0.01, which is consistent with a circular orbit. 
The error is determined from 100 Monte Carlo simulations. The 3a upper bound on the 
eccentricity is e < 0.03. We have also attempted varying the third light contribution to the 
system, flnding L3 = 0.06 ± 0.14 which is consistent with zero, though poorly constrained. 

Combining the parameters determined from the RV curves with the parameters deter- 
mined from the light curve, we flnd the following masses and radii for the stars: Mp = 
1.58 ± O.O4M0, Rp = 1.39 ± O.OTi?©, Ms = 1.58 ± O.O4M0, and Rs = 1.38 ± 0.07i?Q. The 
errors in the masses are dominated by the uncertainties in the velocity semi-amplitudes, 
while the errors in the radii are dominated by the large uncertainty in Rs/Rp. 



In flgure we plot the age-metallicity contours for the eclipsing binary system. To 
calculate at each age/metallicity point we use the Yale-Yonsei version 2 isochrones (Y2 



Demarque et al, 



20041 ) and find the pair of masses which minimize 



2 ^ {Kp-Kp^,soc{Mp,P,siYii)f {Ks - Ks,isociMs,P, sim)y (a - ajsodMp, Ms, P))' 



{k - kisociMp, Ms)f , (T,,e - Ta,e,isociMp, Ms))' 



+ 



where the "isoc" subscript denotes values calculated from the isochrones using the trial 
masses (Mp and Ms) and the measured orbital period (P) and sini, a = {Rp + Rs)/a, 
k = Rs/Rp, and Tave is the average temperature of the two stars. The contours in figure 1281 
show the 68.3%, 95.4% and 99.7% confidence levels. The point shows the values for the 
cluster using the Y2 isochrones (see Paper I). The binary system is younger and has a lower 
metallicity than the cluster. 

Figure [29] shows the observed masses and radii of the components together with the 
expected relation for the cluster from the Y2 isochrones and the relation for [M/H] = —0.25 
and an age of 100 Myr. The stars appear to have radii that are too small given their masses 
for them to have the metallicity and age of the cluster. 

The individual r magnitudes of the stars are Vp = 15.53 and = 15.55. As we stated 
at the beginning of this subsection, the expected masses for the stars if the system were a 
member of the cluster would be Mp = Mg = 1.18 ± 0.04Mq, with an uncertainty dominated 
by the systematic uncertainty in the M — r relation. The observed masses are ~ 7cr above 
the expected masses, this together with the metallicity discrepancy (see figure 128]) . and the 
~ 3cr discrepancy between the systemic velocity and the mean RV of the cluster stars leads 
us to conclude that the system is most likely not a member of the cluster, and is instead 
located in the background galactic disk. 

Finally, we return to the nature of the non-eclipse variability seen in the light curve. 
The periods and amplitudes of the non-eclipse variations are consistent with one or both of 
the stars being a 7 Doradus pulsator. The spectroscopic temperature is slightly lower than, 
though consistent with, the instability strip for this type of pulsator. 

As seen in figure the eclipsing binary components appear to be relatively young 
(< 700 Myr at 3a, < 200 Myr at la). This is consistent with s ome indications th at the 7 



Doradus phenomenon is restricted to stars with ages < 250 Myr (IKrisciunas Ill998l ). though 



there are potentially a few 7 Doradus variables as old as 1 Gyr that have been discovered 



( jPepper and Burke II2OO6I ). To our knowledge, V1380 is only t he second 7 Doradus variable 



found in an eclipsing binary, the other system being VZ CVn (jibanoglu et al. 1120071 ). These 
systems are particularly interesting since the geometry allows for the direct determination 
of the masses and radii of the stars, and with precise photometry it may be possible to use 
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the eclipses to aid in identifying the pulsation modes (IRiazi and Abedi II2006I ). 



6.5. Detached Eclipsing Binary V1028 

Based on its photometry, the detached eclipsing binary V1028 is a candidate cluster 
member. V1028 matches to 2mass05523843+3223296, it has an orbital period of 5.496 ±0.18 
days, a primary eclipse epoch of HJDq = 2453725.6977, and photometry B — V = 1.591, 
V = 19.428, g-r = 1.61, r-i = 0.93, r = 18.73, J - Kg = 0.695, H - Kg = 0.137, and 
Ks = 15.289. If the system is a member of the cluster, its primary would have a mass of 
0.6 < M < 0.7Mq based on the r magnitude of the system. We will now demonstrate that 
the components of this system likely have masses of ~ 0.6Mq and may very well be members 
of the cluster. 

This system shows out of eclipse variations with a peak-to-peak amplitude of 0.028 mag. 
Since the maximum of this variation does not occur at either of the quadrature points or 
near the eclipses, and since the system is well detached, the variation is most likely due to 
spots on one of the stars rather than due to proximity effects. This confirms that at least 
one of the stars has a rotation period synchronized to the orbital period. The primary and 
secondary eclipses occur 0.5 apart in phase, so it appears that the system is also circularized. 

Before fitting an eclipse model to the light curve, we first remove the out of eclipse 
variation using a sinusoid. Assuming the cluster reddening and assuming all the light comes 
from the primary component, we estimate the temperature of the primary is ~ 5000 K from 
the J — Ks and H — colors, f or which the r quadratic limb darkening coefficients are 



a = 0.538, b = 0.199 (IClaret Il2004f ). Using JKTEBOP, we fit the light curve assuming zero 
eccentricity and find the parameters listed in table [91 As before we estimate the parameter 
errors using 100 Monte Carlo simulations taking the error at each observation to be 8.0 
mmag. In this case we refine the period and primary eclipse epoch using JKTEBOP. We 
find that it is necessary to include a significant amount of third light (L3 = 0.266 ± 0.006) 
in the fit, though it is unclear if this is an artifact of not correctly modelling the spot 
through the eclipses. The resulting fractional radii for the primary and secondary stars are 
Rp/a = 0.0439 ± 0.0003 and Rs/a = 0.0416 ± 0.0004. Figure [30] shows the best fit model to 
the light curve. 

To estimate the masses of the component stars we use Kepler's third law and the mass- 
radius relation for low-mass stars. Assuming a power-law mass-radius relation of the form 



- 26 - 



M = Mq{R/ Rq) , the mass of each component is given by: 

b/{3-b) 

Mp,, = Mo(^) |74.471^(— -) I (-) +(-) 1 1 . (10) 




For the mass- radius relation we use the double-hned detached echpsing binary data given in 



Ribas I (120061 ) to find (M/Mq) = 1M{R/ Rq^^ "^^ 



Using eq. [TO]we find component masses of Mp ~ 0.61Mq, and Ms ~ O.57M0. The mass 
of the primary is thus close to what we would expect based on the r magnitude of the system 
if the stars were members of the cluster. The individual magnitudes of the two component 
stars are rp ~ 19.5 and ~ 20.1, with a third component having ~ 20.4. The magnitude 
difference of 0.59 ± 0.03 between the primary and secondary stars is quite a bit larger than 
the value of 0.3 that one would expect from the estimated masses (see the mass-r magnitude 
relation in Paper 1). However, given that we do not properly model the spots on the stars 
through the eclipses, and that we do not have radial velocity data to confirm the masses, 
the discrepancy may simply be an artifact of our modelling. 

We have shown that the stars likely have masses of ~ O.6M0 and may very well be mem- 
bers of the cluster. The system is therefore interesting for follow-up as only ~ 6 double-lined 
detached eclipsing binary systems with main sequence stars smaller than the components in 
this system a re known, and the models for these stars are known to be in error by 5 — 20% 



(iRibas I I2OO6I ) . If the system proves to be a member of the cluster, then the additional 
constraints on its age, metallicity and apparent distance modulus would make the system 
particularly useful in testing low mass stellar models. The downside though to any follow-up 
is that with V ~ 19.4, a significant amount of large telescope time would be needed to obtain 
its radial velocity curve. Note that the expected semi-amplitude of the radial velocity curve 
is ~ 60 km/s. 



7. Conclusion 

In this paper we have identified 1430 new variable stars in the field of M37. This 
includes 20 new eclipsing binaries and 31 new short-period pulsating stars. Approximately 
500 of the variables are F-M main sequence cluster members; for these stars the variability 
is due to rotation coupled with significant surface brightness inhomogeneities. This data set 
is significantly larger than any other set of stellar rotation periods for a cluster that is > 500 
Myr old and will thus provide a unique window on the rotation of lower main sequence stars. 
We argue that many of the other variable stars are young, rapidly rotating main sequence 
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stars in the galactic disk by showing that the amphtude of these variables is anti-correlated 
with the Rossby number. 

We have used 10 W UMa eclipsing binaries and 4 fundamental mode 5-Scutis to inves- 
tigate the relation between distance and extinction along the line of sight, and have shown 
that the extinction apparently cuts off at nearly twice the scale height of the disk. 

Finally, we have analyzed a few particularly interesting variable stars including a previ- 
ously identified multi-mode pulsator which we argue to be a hybrid 7-Doradus/5-Scuti, two 
possible quiescent cataclysmic variables, a DEB with at least one 7-Doradus pulsating com- 
ponent, and a low mass DEB that is a possible cluster member. We have combined radial 
velocity and light curves to obtain a physical solution for the DEB+7-Doradus system. 

This research has made use of the SIMBAD database, operated at CDS, Strasbourg, 
France. We are grateful to C. Alcock for providing partial support for this project through 
his NSF grant (AST-0501681). Funding for M. Holman came from NASA Origins grant 
NNG06GH69G. We would like to thank G. Torres for help with running the eclipsing binary 
analysis software and for helpful comments, G. Fiiresz and A. Szentgyorgyi for help in 
preparing the Hectochelle observations, S. Meibom for help with the stellar classification, 
A. Schwarzenberg-Czerny for discussion of the period finding techniques, J. Kaluzny for 
pointing out the CV interpretation of VI 148 to us and the staff of the MMT, without 
whom this work would not have been possible. We would also like to thank the MMT 
TAG for awarding us a significant amount of telescope time for this project, and the referee, 
D. Weldrake, for thoughtful comments that improved the quality of this paper. 



A. The VARTOOLS Program 

In conducting the research described in this paper we have developed software to perform 
a number of common analysis routines on a set of light curves. The program is called 
Vartools and we release it to the communitj]^. The program reads in a light curve, or a 
set of light curves, in ascii format and performs a series of commands on them; it is designed 
to run the light curves through a pipeline of routines. 

At the time of writing the commands that can be run include the period finding routines 
discussed in §3, routines to fit a Fourier series to light curves, routines to linearly decorrelate 
light curves against generic signals, and implementations of the SYSREM and TEA de- 



^http: / /www. cfa.harvard.edu/~jhartman/vartools 
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trending algorithms mentioned in §3.1. This is not an exhaustive hst; a full description of 
the program can be found on the website. 
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Fig. 1. — RMS of the light curves as a function of r magnitude. The point-to-point precision 
of the Hght curves approaches 1 mmag at the bright end. 
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Fig. 2. — Histogram of peak frequencies in the Lomb-Scargle periodograms of all light curves. 
The top panel shows the histogram before removing bad frequencies (seen as spikes in the 
periodogram) , the bottom shows the histogram after removing them. 



-35- 




Fig. 3. — Histogram of the logarithm of the formal false alarm probability values from the 
Lomb-Scargle period finding algorithm. Example phased light curves at a range of FAP 
values are shown. The units on the inset plots arc phase for the x-axis and r magnitude 
for the y-axis. Wc select stars with FAP < —150 and P > 0.1 days or FAP < —10 and 
P < 0.1 days as variables. The dotted line shows the selection for the longer period variables. 
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Fig. 4. — Period vs. the negative logarithm of the formal false alarm probability values from 
the Lomb-Scargie period finding algorithm. We only plot points the pass the removal of 
bad frequencies described in the text and that have FAP < —2. The dotted hue shows the 
selection. Dark points show selected variable candidates, light points show rejected stars. 
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Fig. 5. — Same as figure [2], this time shown for AoV. 
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Fig. 6. — Same as figure [3l here we plot the histogram of AoV. The selection criterion is 
AoV > 3.75. 
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Fig. 7. — Same as figure [2], this time shown for BLS. 
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Fig. 8. — Same as figure [3], this time shown for BLS. The selection criterion is SN > 12.0. 
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Fig. 9. — The uncertainty on the variable star periods as a function of period. The sohd hues 
shows the expected relation for a signal to noise of 10 and a residual correlation time-scale 
of 30 minutes. 
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Fig. 10. — Examples of candidate variable star light curves which were rejected by eye. The 
listed period is the value returned by LS and is in days. The top candidate may very well 
be real, we erred on the side of caution, however, and rejected it since it is bright enough 
for the light curve to be dominated by red noise, shows signficant scatter when phased and 
lies close to the edge of the chip where the image subtraction tended to be less reliable. The 
next two candidates are both located on the same chip. The apparent variation is due to 
the same five nights being brighter than average for both stars. This systematic variation 
was seen in the light curves of several stars on this chip and is due to the appearance of bad 
columns and rows on the chip that did not occur on other nights. The second candidate 
from the bottom may be real but was rejected due to its proximity to the edge of the chip. 
The bottom candidate was rejected because it lies in the wings of a much brighter saturated 
star. 



-43- 



V3 P=0,43a5ESK 
L&.6 



V4 P=0,5MST3Sfl 



■V5 P=0,ffi'e794efi 



V9 p=o,ioee3M3 

15,7 



'-I 16.63 




Fig. 11. — Light curves of 15 previously identified variabl es tfiat we naatcfi to. For vari- 
ables V3-y 7 we ov erplot the Rc published photometry from iKiss et al. I (l200ll ) in gray. The 
Kiss et al. I (120011 ) photometry has been shifted by an offset to match up with our data. 
Note that V4 exhibits a pronounced O'Connell effect and also shows evidence for variations 
of ~ 0.01 mag between cycles. 
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Fig. 12. — Location of all variables identified by this survey on gr and gi CMDs. The dark 
points are the variables, the gray-scale is used to show all the point sources. The variables 
appear to lie preferentially along the cluster main sequence. The cut-off at the bright end is 
due to saturation of our primary time-scries data. Note that the formal photometric error 
on the colors goes from ~ 0.015 mag at the bright end of each plot to ~ 0.03 mag at the 
faint end. 
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Fig. 13. — Location of visually selected short-period pulsating variables (filled points) and 
eclipsing binaries (open squares) on gr and gi CMDs. Labels indicate some of the variables 
discussed in the text. Note that we do not use a symbol to plot V898 as it is neither an 
eclipsing binary nor a pulsating star, the label however shows its location on the CMDs. 
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Fig. 14. — Phased light curves for 31 visually selected short-period pulsating variables. Note 
that several of these stars, including V72, V108, V1244 and V1368, are multi-periodic. 
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Fig. 14. — Continued. 
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Fig. 15. — Phased light curves for 20 visually selected eclipsing binaries. 
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Fig. 15. — Continued. 
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r [Mag] 

Fig. 16. — Period-r relation for (left) all variables stars, (center) photometrically selected 
candidate cluster members and (right) non-cluster members. Note that there is a strong 
correlation between period and magnitude for cluster members while for non-members there 
is no clear correlation. 
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Fig. 17. — Ay vs distance for 10 W UMa eclipsing binary systems (open squares), 4 S- 
Scuti stars (filled p oints) and the open cluster (star). The line shows the Besangon model 
jRobin et al. Iboosh : or the galactic latitude/longitude of the field, assuming an interstellar 
extinction of 0.7 mag/kpc. 
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Fig. 18. — Phased light curves for a random sample of 16 variable stars that are not members 
of the cluster, eclipsing binaries, or short-period pulsators. Many of these variables are likely 
young, heavily spotted, rapidly rotating stars. 
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Fig. 19. — Period vs. V — Ic for variable stars that are not members of the cluster, not 
eclipsing binaries, and not 5-Scuti or RR Lyrae type pulsators. The left panel shows the 
period on a linear scale, the right panel on a logarithmic scale. 
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Fig. 20. — B — VvsV — Ic color-color plot showing all stars (gray points), field variables that 
are not eclipsing binaries nor obvious pulsators (dark points), and the fiducial de-reddened 
relation for the cluster (sohd fine). The arrow shows the direction of the reddening vector. 
Using this diagram we estimate the reddening to all the stars with B — V< 1.2, V — Ic< 1-8. 
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Fig. 21. — The r-band amplitude (A^) versus the Rossby number (i?o) for field variables not 
classified as eclipsing binaries or pulsating stars (see text for a description of the selection). 
The filled points show the measured values, the open circles show the minimum amplitude 
that each variable could have had and still have been detected. The anti-correlation between 
Ro and is consistent with the hypothesis that these variables are heavily spotted, rotating 
stars. 
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Fig. 22. — The DFT power spectrum of the hght curve for V17. The top panel shows the 
window power spectrum, the second panel shows the raw DFT power spectrum, the third 
panel shows the power spectrum after implementing the CLEAN deconvolution algorithm, 
the bottom panel shows the CLEAN power spectrum zoomed in on the low frequency region. 
The bottom three plots have been scaled by a factor of 10^. The CLEAN power spectrum 
reveals five frequencies at 22.796 days'S 1.418 day-\ 2.263 day-\ 2.054 day'S and 0.823 
day~^. 
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Fig. 23. — Light curve of the variable V898. The hght curve is non-periodic, showing 
continuous variabiUty with repeated outbursts having amphtudes of 0.1 — 0.2 mag, durations 
of 0.5 days and a time between outbursts of 1-2 days. 
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Fig. 24. — The DFT power spectrum of the hght curve for the detached echpsing binary 
system V1380. The top panel shows the window power spectrum, the middle panel shows the 
raw DFT power spectrum, the bottom panel shows the power spectrum after 10 iterations 
of the CLEAN deconvolution algorithm. The bottom two plots have been scaled by a factor 
of 10^. The CLEAN power spectrum reveals two strong signals, one at a frequency of 
L063 ± 0.013 days^^, the other at a frequency of 0.852 ± 0.013 days^^. These correspond to 
periods of 0.941 ± 0.011 days, and 1.174 ± 0.013 days, the same periods found by L-S. 
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Fig. 25. — Light curve of the detached echpsing binary system V1380. Top: phased at the 
orbital period of 2.1916 days. Second from top: out of echpse points phased at 0.941 days. 
Second from bottom: out of echpse points, after subtracting a sinusoid fit to the second from 
top plot, phased at 1.174 days. Bottom: phased at the orbital period after removing the 
signals in the middle two panels. 
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Fig. 26. — Radial velocity curves for the primary and secondary components of the double- 
lined eclipsing binary system V1380. The solid line shows a zero-eccentricity fit to the 
observations. We assume an error of 2.2 km/s for each observation. The stars have equal 
mass to within 1%. 
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Fig. 27. — Best fit detached eclipsing binary model to the light curve of V1380. The top 
panel shows the full phased light curve, the bottom panels are zoomed in on the echpses. 
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Fig. 28. — Age-metallicity contours for the eclipsing binary system V1380 from comparing 
the Y2 isochrones with the observed masses/radii/temperatures of the components. The 
contours show the 68.3%, 95.4% and 99.7% confidence levels. The point shows the values 
for the cluster. The binary system appears to have a lower metallicity and younger age than 
the cluster. 
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Fig. 29. — The masses and radii of the two components of V1380 together with the expected 
relations from the Y2 isochrones for the cluster metallicity of [M/H] = 0.045 and age of 580 
Myr (dotted line), and for a lower metallicity of [M/H] = —0.25 and age of 100 Myr (solid 
line). The stars appear to have radii that are too small, given their masses, for them to have 
the metallicity and age of the cluster. 



-64- 




Fig. 30. — Best fit detached ecfipsing binary model to the fight curve of V1028. The top 
panel shows the full phased light curve, the bottom panels are zoomed in on the eclipses. In 
this plot we adopt the zero-point of phase to be the time of primary eclipse. 



Table 1. Variable Star Catalog 



ID RA (J2000) DEC {J2000) Chip Period (days) crp{days) < r > (mag)''' iJA/Sr (mag)'' s (mag)'= r (mag)"^ i(mag)'= S(mag)'^ V (mag)'' Type"-' Selection' 



1 


05:52:20.35 


+32:33:20.4 


















14.015 


13.406 


10 


000 


2 


05:52:16.54 


+32:28:15.7 


















15.341 


14.653 


10 


000 


3 


05:52:33.01 


+32:32:42.0 


14 


0.42248314 


0.00001000 


15.7511 


0.1057 


16.350 


15.789 


15.460 


16.779 


16.051 


10 


111 


4 


05:52:53.26 


+32:33:01.5 


5 


0.55819369 


0.00002600 


15.6407 


0.0815 


16.396 


15.634 


15.157 


16.886 


15.948 


10 


111 


5 


05:53:00.64 


+32:24:50.8 


8 


0.27878242 


0.00000400 


15.9035 


0.1140 


16.374 


15.915 


15.433 


16.769 


16.184 


01 


111 


6 


05:51:50.53 


+32:32:34.8 


23 


0.10983967 


0.00000100 


15.9552 


0.1263 


16.564 


15.973 


15.800 






01 


111 


7 


05:52:39.09 


+32:36:31.2 


13 


0.35773482 


0.00000400 


17.3113 


0.1562 


18.290 


17.252 


17.195 


18.705 


17.890 


10 


111 


8 


05:52:34.32 


+32:32:18.7 




0.11950000 














12.925 


12.498 


01 


000 


9 


05:52:14.91 


+32:24:41.0 




0.09085950 














13.668 


13.298 


01 


000 


10 


05:52:00.49 


+32:36:48.2 




0.07843140 














15.293 


14.762 


01 


000 


11 


05:52:03.38 


+32:35:13.7 


22 


0.11867118 


0.00000200 


15.2542 


0.0296 


15.669 


15.245 


15.059 


15.908 


15.456 


01 


110 



^Thc flux averaged r magnitude of the light curve. 
^Thc root-mean-square of the light curve. 

'^Thc magnitude of the source in the photometric catalog (Paper I) 
Value from lKalirai et al. I i200lf) 

two-bit flag. The first bit denotes whether or not the star is an eclipsing binary. The second denotes whether or not it shows pulsations. 

three-bit flag. The first bit denotes if the variable was selected by LS, the second if it was selected by AoV and the third if it was selected by BLS. 



Table 2. Summary of Variable Star Catalog 



Number 



Total 1454 

Previously Known 24 

Total Detected in Survey 1445 

Recovered Previously Known 15 

New Discoveries 1430 

Eclipsing Binaries 31 

Pulsating Variables 44 

With BVgri Photometry 1378 

Near gri main sequence 753 

Near BVgri main sequence 527 



Table 3. Catalog of Suspected Variable Stars 



ID RA (J2000) DEC {J2000) Chip Period (days) CTp (days) < r > (mag)" RMSr (mag)'' g (mag)= r (mag)= i (mag)"^ B (mag)"* V (mag)"* Type"-' Selection'' 



1 


05 


51 


20 


85 


+32 


34 


50 


5 


31 


7 


17495773 





02262700 


18 


2823 





0189 


18 


217 






20 


295 


18 


923 


00 


110 


2 


05 


51 


29 


80 


+32 


35 


42 





31 










14 


4660 





1214 


14 


466 






14 


699 


14 


355 


00 


001 


3 


05 


51 


33 


12 


+32 


30 


33 


8 


33 


6 


50096134 





01918200 


19 


4794 





0157 


19 


479 


18 


653 


21 


571 


20 


077 


00 


100 


4 


05 


51 


37 


12 


+32 


40 


23 


4 


29 


4 


55114850 





00475400 


15 


8300 





0037 


15 


829 


15 


537 


16 


886 


16 


071 


00 


100 


5 


05 


51 


37 


70 


+32 


40 


40 


7 


29 





85990601 





00037500 


21 


4176 





0440 


21 


383 














00 


100 


6 


05 


51 


42 


45 


+32 


35 


40 


4 


31 


15 


66710882 





11219700 


18 


9369 





0084 


18 


939 


18 


418 


21 


483 


20 


447 


00 


100 


7 


05 


51 


45 


62 


+32 


27 


16 


7 


34 


11 


85243048 





08133600 


19 


5583 





0205 


19 


564 


18 


635 


21 


824 


20 


221 


00 


100 


8 


05 


51 


50 


79 


+32 


34 


01 


8 


23 





67217666 





00018800 


20 


2537 





0371 


20 


282 


19 


141 


22 


477 


21 


063 


00 


110 


9 


05 


51 


51 


70 


+32 


38 


43 





21 





29449917 





00007200 


23 


5355 





2779 


23 


506 


22 


150 










10 


010 


10 


05 


51 


56 


52 


+32 


43 


44 


1 


19 


11 


73852830 





03931500 


18 


3177 





0081 


18 


318 


19 


454 


20 


049 


18 


656 


00 


100 


11 


05 


51 


56 


83 


+32 


37 


02 


5 


22 





60196052 





00019900 


20 


8365 





0395 


20 


824 


19 


529 


22 


803 


21 


452 


00 


100 


12 


05 


51 


57 


09 


+32 


33 


08 


4 


23 





89557446 





00023500 


22 


4039 





1092 


22 


387 


20 


483 


24 


875 


23 


399 


00 


100 


13 


05 


51 


59 


78 


+32 


23 


33 


2 


26 





93428743 





00049300 


19 


2657 





0101 


19 


266 














00 


100 


14 


05 


52 


00 


97 


+32 


28 


10 





25 


1 


37550745 





00056300 


17 


0477 





0123 


17 


049 


16 


885 










00 


100 


15 


05 


52 


03 


88 


+32 


39 


10 


1 


21 


2 


85961673 





00233400 


19 


2151 





0211 


19 


207 


18 


404 


21 


225 


19 


879 


00 


100 


16 


05 


52 


09 


58 


+32 


26 


58 


4 


25 


18 


41735750 


3 


59827300 


22 


7522 





1035 


22 


770 


20 


951 


25 


276 


23 


538 


00 


010 


17 


05 


52 


15 


33 


+32 


35 


44 


7 


22 





74276537 





00021900 


21 


0196 





0436 


21 


Oil 


19 


582 










00 


100 


18 


05 


52 


15 


42 


+32 


3D 


46 


7 


24 





90985227 





00054200 


21 


9206 





1236 


21 


956 


20 


277 










00 


100 


19 


05 


52 


16 


06 


+32 


32 


44 


6 


23 


5 


25393598 





01656200 


21 


0112 





0550 


21 


018 


20 


086 










00 


100 


20 


05 


52 


17 


60 


+32 


30 


12 


6 


24 


8 


46159122 





01841700 


17 


9894 





0125 


17 


983 


17 


305 


19 


724 


18 


465 


00 


100 


21 


05 


52 


30 


27 


+32 


3D 


49 


8 


15 


1 


42595003 





00214900 


18 


8178 





0092 


18 


822 


18 


356 


19 


981 


19 


118 


00 


100 


22 


05 


52 


34 


04 


+32 


29 


22 


6 


15 


8 


58126897 





01497500 


17 


9452 





0057 


17 


949 


17 


591 


19 


033 


18 


214 


00 


o^io 


23 


05 


52 


34 


50 


+32 


4D 


23 





11 


2 


09687226 





00201000 


20 


2432 





0209 


20 


218 














00 


^100 


24 


05 


52 


40 


10 


+32 


36 


10 


4 


13 


11 


89564092 





07649400 


18 


9641 





0087 


18 


970 


18 


218 


21 


161 


19 


489 


00 


100 


25 


05 


52 


40 


60 


+32 


31 


50 


1 


14 


4 


47741049 





00336800 


19 


0914 





0349 


19 


088 


18 


364 










00 


110 


26 


05 


52 


50 


60 


+32 


34 


48 


9 


4 


5 


20926065 





19755000 


15 


2892 





0151 


15 


290 


14 


945 










00 


110 


27 


05 


52 


54 


61 


+32 


39 


28 


6 


3 


17 


23385097 





08937600 


15 


8058 





0052 


15 


808 


15 


450 










00 


100 


28 


05 


53 


09 


81 


+32 


21 


34 


5 


9 


2 


59599006 





00113400 


15 


6733 





0148 


15 


655 






16 


310 


15 


699 


00 


111 


29 


05 


53 


11 


99 


+32 


23 


53 


9 


8 





58837547 





00003100 


18 


2538 





0971 


18 


298 


17 


203 










10 


111 



^Thc flux averaged r magnitude of the light curve. 
'^Thc root-mean-square of the light curve. 

'^The magnitude of the source in the photometric catalog {Paper I). Note that when the i magnitude is missing no color term was applied in transforming from the instrumental r and i magnitudes to 
the Sloan 2.5 m system. 

"^Value from lKalirai et al. I i200l|) 

'^A two-bit flag. The first bit denotes whether or not the star is an eclipsing binary. The second denotes whether or not it shows pulsations. 

three-bit flag. The first bit denotes if the variable was selected by LS, the second if it was selected by AoV and the third if it was selected by BLS. 
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Table 4. Extinction and Distances for Fundamental Mode (5-Scuti Stars 



ID 


V 


B-V 


V-Ic 


Period (days) 


Av 


Distance (kpc) 


V5 


16.184 


0.585 


1.194 


0.278782 


1.56 


7.90 


V6 


16.177 


0.691 


0.806 


0.109840 


0.95 


5.22 


Vll 


15.456 


0.452 


0.826 


0.118671 


0.98 


3.92 


V840 


20.695 


1.711 


1.377 


0.145197 


2.04 


31.21 
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Table 5. Extinction and Distances for W UMa Eclipsing Binaries 



ID 


V 


B-V 


V-Ic 


Period (days) 


Av 


Distance (kpc) 


V3 


16.051 


0.728 


1.026 


0.422483 


1.18 


2.78 


V4 


15.948 


0.938 


1.234 


0.55819 


1.27 


2.04 


V7 


17.890 


0.815 


1.117 


0.357735 


1.23 


5.10 


V20 


17.628 


0.892 


1.008 


0.289456 


0.45 


3.07 


V24 


20.506 


1.275 


1.518 


0.252674 


1.23 


5.84 


V855 


18.156 


1.227 


1.617 


0.564829 


1.86 


3.56 


V1160 


20.234 


1.134 


1.435 


0.242092 


1.42 


6.88 


V1181 


21.847 


1.402 


1.871 


0.25817 


2.32 


10.30 


V1194 


19.030 


1.136 


1.267 


0.269735 


0.56 


3.58 


V1447 


19.930 


1.066 


1.535 


0.296877 


2.13 


8.52 
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Table 6. Radial velocity and light ratio measurements from the spectra of V1380. 



HJD Primary Velocity [km/s] Secondary Velocity [km/s] Light Ratio [L2/ii] 



2454154.6776 -104.36 115.64 0.991 

2454165.7104 -108.62 119.24 0.955 

2454170.6082 -31.08 40.47 0.995 

2454171.6861 49.34 -38.14 1.025 
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Table 7. Eclipse time measurements for V1380. The last measurement given is inferred 

from a fit to the radial velocity data. 



Eclipse Number 


Observed Minimum [HJD] 


Error [HJD] 


0.5 


2453725.6289 


0.0008 


1.0 


2453726.7223 


0.0004 


1.5 


2453727.8188 


0.0003 


2.0 


2453728.9160 


0.0004 


7.5 


2453740.9756 


0.0003 


195 


2454152.0803 


0.0037 
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Table 8. Parameters for Eclipsing Binary V1380. 



Parameter 


Value 


Error 


■Porb 


2.19258 days 


0.00004 days 


HJDo 


2453724.5306 


0.0004 


Ki 


118.3 km/s 


1.3 km/s 


K2 


118.3 km/s 


1.3 km/s 


7 


5.3 km/s 


0.7 km/s 


Js/Jp 


1.02 


0.02 


(Rp + Rs)/a 


0.266 


0.002 


Rs/Rp 


0.99 


0.1 


i 


80.00° 


0.09° 


a 


10.41 Rq 


0.23 Rq 


Ml 


1.58 Mq 


0.04 Mq 


M2 


1.58 Mq 


0.04 Mq 


Ri 


1.39 Rq 


0.07 Rq 


R2 


1.38 Rq 


0.07 Rq 
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Table 9. Light Curve Parameters for V1028. 



Parameter 


Value 


Error 


Porh ~ Prot 


5.496 days 


0.18 days 


HJDo 


2453725.6989 


0.0001 


Js/Jp 


0.649 


0.004 


{Rp + Rs)/a 


0.0855 


0.0003 


Rs/Rp 


0.95 


0.01 


i 


89.73° 


0.05° 




0.266 


0.006 



